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1.  Introduction 

Recently  I  have  spent  some  effort  on  a  grant  supported  by  the  U.S.  Army 

Materials  and  Mechanics  Research  Center  in  studying  numerically  and  analytically  the 

diffusion  of  agent  from  a  sessile  agent  droplet  (a  spherical  cap)  through  the 

1-4 

(polymer)  membrane  on  which  it  rests.  The  agent  droplet  can  be  neat  or  polymer 

thickened. Dr.  Nathaniel  Schneider  of  Natick  R.  D.  and  E.  Laboratories  Who 

has  asked  us  to  review  these  matters  is  primarily  interested  only  in  the  neat  agent 

problem.  This  diffusion  problem  arises  in  assessing  the  definite  amount  of  agent, 

Q^,  coming  through  a  test  membrane  within  a  definite  time  period,  t^.  Its 

relevance  to  monitoring  and  dosimetry  is  thus  clear  and  needs  no  further  comment. 

The  numerical  solutions  make  clear  both  quantitatively  and  quantitatively  (for 

suitably  chosen  parameters)  the  overall  expected  behavior  and  the  results  are  in 

overall  agreement  with  observations.  While  for  the  polymer  thickened  agent  droplet 

problem  there  is  a  more  extensive  region  of  scaling  with  regard  to  initial  droplet 

radius,  R^,  and  membrane  thickness,  L,  this  is  not  the  case  for  the  neat 
1  1-3 

droplet.  Specifically  we  showed  that  simple  scaling  relations  apply  only 
for  a  rather  limited  parameter  ranges:  a)  The  more  ubiquitous  (under  field 
conditions)  case  where  the  agent  is  deposited  as  an  aerosol  of  tiny  droplets  tihich 
spread  on  the  membrane,  for  which  R^/L  «  1,  the  small  droplet  case.  This  case  is 
essentially  coiiq>letely  analyzed  in  references  (1)  and  (4)  and  a  quasi  steady  state 
analytical  theory  is  given  in  reference  (2)  for  the  case  that  the  agent  solubility 
is  sufficiently  large.  I  am  given  to  understand  by  Dr.  Schneider  that  the  scaling 
and  applicability  of  this  case  is  sufficiently  well  understood  that  no  further 
comment  is  necessary  here,  (b)  The  more  special  case,  we  call  the  large  droplet 
case,  arises  when  R  /L  »  1.  Note  that  the  condition  used  for  approximate 
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solution  of  this  problem  that  the  droplet  is  smaller  than  L  cannot  be  correct  for 
all  time.^  The  radius  of  the  sessile  droplet  decreases  continuously  with  time  and 
will  ultimately  become  equat  to  or  smaller  than  L  (as  the  agent  diffuses  through  the 
membrane)  when  the  special  scaling  for  the  limiting  ease  must  fail.  In  the  next 
section  we  review  the  discussion  of  this  special  limiting  case.  It  has  been 
originally  treated  by  us  in  reference  (1)  and  (3).  We  shall  treat  it  here  more 
extensively  in  the  next  section  showing  that  it  arises  from  a  singular  perturbation 
of  the  original  boundary  value  problem.  In  the  last  section  we  deal  with  the 
extension  of  this  problem  when  evaporation  from  the  sessile  droplet  surface  is  not 
negligible. 

2.  The  large  droplet  case  -  negligible  surface  evaporation  from  the  droplet. 

We  consider  the  moving  boundary  diffusion  problem  initially  described  in 

references  (1)  -  (4).  We  use  a  slightly  modified  notation.  All  physical  entities 

with  dimensions,  other  than  certain  scaling  variables  such  as  R  ,  L,  etc.  will  be 

0 

written  with  a  carat  over  the  letter  denoting  the  entity.  We  will  subsequently 

introduce  suitably  scaled,  dimensionless  variables  denoted  by  letters  without  carats. 

The  concentration  c  of  agent  in  the  membrane  at  location  r,  z  at  time  t 

1-4 

satisfies  the  diffusion  boundary  value  problem 

Oc/at)  -  D  {(a2c/ar2)  +  l/r  Oc/dr)  +  ibK/bz^)}  (1) 
inO<z<L  ,  r>0,t>0 

with  initial  and  boundary  conditions 

c(r,  z,  t  *  0)  a  0,  (2) 

A  A  A  ^  ^ 

c(r,  0,  t)  =  c  ,  r  <  R(t)  (3) 

o  ~ 

with  c^  the  solubility  of  agent,  i.e.  the  equilibriiim  concentration  of  agent  in  a 
membrane  immersed  in  neat  liquid  agent  at  ambient  fixed  temperature  and  pressure  and 

A  A 

R(t)  the  time  dependent  sessile  droplet  radius.  The  sessile  droplet  is  a  spherical 
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cap  with  fixed  contact  angle  0  and  voliune  1/3  f(0)R^(t), 

_2 

f(0)  =  irsin0(2+cos0)/(l+cos0)  Furthermore, 

(3c/az)A_Q  =  0  ,  r  >  R(t)  (4) 

and 

c(r,  L.  t)  =  0.  (5) 

A  A  A 

The  moving  boundary  of  the  spherical  cap  r  =  R(t)  is  determined  by  the  mass  balance 

A  ^ 

R{  t) 

pf(0)R2  (dR/dt)  =  2irDjdr  r  (Bc/dz)^^Q  ;  R(t=0)  =  Rq  ,  (6) 

0 

with  p  the  density  of  neat  agent.  We  now  introduce  dimensionless  variables: 

A  A 

r  =  r/R  ,  z  =  z/L, 
o 

t  =  Dt/L^  ,  R(t)  =  R(t)/R  , 

o 

c(r,z,t)  =  c(r,S,t)/c  , 

o 

B  =  2irc„/pf(0)  ,  (7) 

o 


and  a  dimensionless  parameter  of  smallness, 

c  *  L/R  «  1, 
o 

into  this  boundary  value  problem. 

One  finds 

Oc/at)  =  (a^c/az^)  +  [<a^c/ar^)  +  i/r  ac/ar)]  (8) 

inO<  z<  1,  r>0,  t>0; 

c(r,  z,  0)  s  0  ,  (9) 

c(r,0,t)  »  1  for  r  <  R(t),  (10) 

(Bc/az)  -  .  0  for  r  >  R(t)  ;  (11) 

ZsU 

c(r,l,t)  -  0  (12) 


and 
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R(t) 

r2  (dR/dt)  =  c  S  J  dr  r  (ac/az)z=0  ♦ 

0 

R(0)  =  1. 

The  natural  way  to  seek  an  approximate  solution  of  (8)  -  (13)  is  to  consider  the 
singular  perturbation  in  which  the  zeroth  order  solution  for  the  concentration 
c^  s  c^(z,t)  satisfies  (8)  with  neglect  of  the  term  in  the  square  bracket  on  the 
r.h.s.  of  (8),  i.e. 

c  =  c^(z.t)  +  0(c)  (14) 

0 

R(t)  =  R^(t)  +  cR^(t)  +  0(c) 
where 

ac  /at  =  a^c  /az^  ,  o<z<i  ,t>0; 
o  o 

c  (z.O)  »  0  ,  (15) 

o 

c^(0,t)  =  1  , 

c  (l,t)  a  0; 
o 

and  using  (14)  in  (13) 

dR  /dt  =  0  ,  R  =1  ,  (16) 

o  o 

®0=1 

(dRi/dt)  a  flj  dr  r  (aco/az)o  =  B/2  (aco/dz)o  ,  (17) 

0 

Rl(0)  =  0 

The  zeroth  order  solution  of  the  boundary  value  problem  (cf.  (14)  and  (15))  has  been 
given  by  us  in  reference  (1)  and  also  in  reference  (3)  where  it  is  given  by  the 
special  case  described  by  Eq.  [7]  on  p.  284  which  in  our  notation  reads: 

Co(z,t)  =  r  {erfc  (2n  +  -  erfc  I(2(n+l)-z/2 /t)  1 )  (18) 

n*l 

t 

R(t)  a  1  +  (c6/2)  /  dt’  ((dCo(z,t’)/dz)] 

0 


ZaO  • 


(19) 
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The  integrals  given  above  are  tabulated  in  standard  tables. 

To  obtain  a  perturbation  solution  kdiich  holds  for  longer  times  one  must 
introduce  "stretched”  coordinates,  and  find  the  "inner"  and  "outer"  solution  of  this 
singular  perturbation.  We  shall  not  carry  out  this  rather  dull  but  lengthy  standard 
problem  in  applied  mathematics. 

Instead  we  point  out  that  the  zeroth  order  solution  provides  a  limiting  case, 
valid  for  very  short  times,  when  we  further  can  replace  the  infinite  series  in  (18) 
by  its  leading  term  for  which 

(aco/az)2=o~  .  (20) 

Substituting  this  into  (17)  yields 

R^(t)  =  -flt(t/ir)^^^  +  {2(t/v)^''^e"^^^  - 

2  erfc  (l//t)}]  (21) 


or 

R(t)  ~  1-  cBl(t/ir)^^^  +  {2(t/ir)^^V^^^  -  2  erfc(l/ /t) }  ] . 

In  zeroth  approximation  the  flux  through  themembrane  at  z  =  1,  (which  in  this 
approximation  is  exactly  the  same  as  the  area  averaged  flux  over  the  membrane)  can 
be  obtained  from  (18).  We  shall  again  only  retain  the  leading  order  term  in  the 
summation  (18)  to  calculate  the  zeroth  order  flux,  which  we  write  in  dimensional 


form, 


F(t)  .  -D  Oc  /az)/'  , 

O  ZsL  ^ 

o  i.Nl/2  -lf/4Dt 

a.  2c  (D/*t)  c  (22) 

o 

The  absence  of  R  in  (22)  clearly  shows  that  (22)  can  only  apply  for  short 
o 


enough  times  even  if  c  <  <  1.  The  use  of  (22)  to  calculate  the  reduced, 

^  2 

dimensionless  break  through  time,  t.  D/L  ,  requires  that 

D 

A  2 
t.  D/L  <  T, 

D  L 


(23a) 


I 
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where  t  is  the  time  needed  for  the  droplet  radius  to  decrease  by  L,  i.e. 

La 

R(t  )  =  1-c .  In  zeroth  order  approximation  the  time  t  is  defined  by  the 

IL  La 

solution  of  the  transcendental  equation 

S  -2  erfc(l/  ^)}] 

=  1.  (23b) 

The  special  scaling  in  (22)  also  noted  by  Angelopoulos,  Schneider  and  Neldon 
subsequently  to  us,  in  an  unpublished  report  applies  only  well  for  sufficiently 
short  times,  even  for  Rq/L  >  i* 

3.  The  zeroth  order  solution  of  the  large  droplet  case  for  non-negligible 
evaporation. 

As  shorn  in  reference  (4)  if  evaporation  is  not  negligible  (4)  has  to  be 
replaced  by  (3c/az)A_-  +  h  [c(r,0,t)-Ul  =  0  ,  r  >  R(t)  (24) 

and  (b)  by 


A  A 


0 


-R2f(0)Yp 

A 

with  h  the  evaporation  coefficient  characterizing  the  flat  polymer-vapor  interface, 

^  A 

U  the  equilibrium  concentration  of  agent  in  the  vapor  above  the  polymer,  y  the 
(per  unit  surface  area)  evaporation  coefficient  of  agent  from  a  neat  agent  surface 
at  ambient  conditions.  In  the  zeroth  order  approximation  we  neglect  all  diffusion 

A  A  ^ 

for  r  >  R(t)  so  that  the  solution  c^(z,t)  is  unaffected  by  (24)  and  is  still  given 
by  (18).  Denoting  by  y  the  appropriately  scaled,  dimensionless  y  defined  by 
y  *  Ly/D 


(26) 
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we  can  vrrite  the  dimensionless  equation  (25)  as 

R 

r2  (dR/dt)  =  c{Bjdr  r  (dc/Bz)g^o  -yR^}  (27) 

0 

From  (27)  and  (14),  Ro(t)  =  1,  but 
1 

(dR^/dt)  =  8/  dr  r  OCo/az)2^0  -2yRi  (28) 

0 

with  Ri(0)  =  0.  The  solution  of  (28)  is 
t 

Rl(t)  =  Bjdf  e-2Y<t-f)  (aco(r,z,f)/3z)2=0  (29) 

0 

which  can  be  evaluated  easily  for  short  times  using  (20),  the  resulting  integrals 
are  all  tabulated  . 
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